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We present the results of a microscopic calculation of the longitudinal thermal conductivity of quasiparticles, 
k xx , in a 2D d-wave superconductor in the vortex state. Our approach takes into account both impurity scattering 
and a contribution to the thermal transport lifetime due to the scattering of quasiparticles off of vortices. We 
compare the results with the experimental measurements on high-T c cuprates and organic superconductors. 



In the last few years measurements of the ther- 
mal conductivity have become one of the most 
powerful probes of unconventional superconduc- 
tors. Materials with linear nodes in the energy 
gap have been predicted to exhibit a residual lin- 
ear in the temperature, T, term T -> 0, 
which is only weakly sensitive to impurity con- 
centration. jl|] Both properties have been firmly 
established in the high-T c cuprates, and a lin- 
ear term has been resolved in the organic super- 
conductor k-(BEDT-TTF) 2 Cu(NCS) 2 .@ 

Information on near-nodal quasiparticles can 
be obtained from the analysis of transport prop- 
erties in the vortex state. The number of un- 
paired quasiparticles increases in an applied mag- 
netic field, and their existence leads to 
the sublinear, in H, increase of k xx (H) below 
T ~ 1K.|^,^| In contrast, at higher temperatures 
the thermal conductivity of the cuprates has been 
found to decrease with the applied field, and of- 
ten becomes nearly field independent above a few 
Tesla. [Q Recently we have proposed a unified de- 
scription of both features in the cuprates. ]8| In 
these materials the experiments are restricted to 
fields H <C H C 2, while in the organic materials 
the entire H — T range can be probed, making 
the analysis of k xx {H,T) important. 

We consider a clean (I ^> £ , where I is the 
mean free path, and £o is the coherence length) 2- 
dimensional d-wave superconductor with a cylin- 
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drical Fermi surface, model the vortex lattice by 
an Abrikosov-likc solution, and use the quasiclas- 
sical method with the approximation scheme of 
Ref.g. For I > {2eH)- 1 / 2 = A the thermal con- 
ductivity is given bym 
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where Nq is the normal state density of states, 
v is the Fermi velocity, P(u> n ) = [1 — 
i v / ¥(2AA/u)V'(u Il )cos 2 2</.]- 1 / 2 , the function 
W(u) — e~ u erfc(— iu), u n = 2iui n Asgn(uj n )/v, 
the renormalized frequency iuj„ = iu> n + ia{iuj n ), 
depends on the self-energy due to impurity scat- 
tering, a, and A is the spatial average of the 
amplitude of the order parameter. The second 
term in the effective scattering rate is due to 
scattering off of the vortex lattice. This term 
has the same symmetry as the superconducting 
gap, and therefore is more important at higher 
temperatures. At low T and H the main ef- 
fect of the magnetic field is to increase the num- 
ber of the excited quasiparticles near the nodes, 
where vortex scattering is weak; consequently the 
thermal conductivity increases. At higher fields 
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the number of quasiparticles is controlled by the 
temperature, the main effect of the vortices is 
to introduce a new scattering mechanism, and 
k xx decreases to a plateau-like feature at fields 
such that v/A > T a(uj = 0), in agree- 
ment with the experimental results on high-T c 
cuprates.[|) The parameters of the theory are the 
dimcnsionless quantity S = AA/v, and the scat- 
tering rate in the normal state, T; the scatter- 
ing phase shift is fixed to either the strong (uni- 
tarity) or the weak (Born) limit. At the mean 
field level S = a[{H c2 - ff)/Hff c2 ] 1/2 where a = 
(A /hv)(hc/2e) 1/2 . For high-T c materials a w 2- 
7T 1 / 2 , § while for k-(BEDT-TTF) 2 Cu(NCS) 2 
using Ao = 2.4meV and v = 5 ■ 10 6 cm/s || we 
find a ~ 1.3T 1 / 2 . 

At low T in this material k xx (H) decreases 
sharply below -ff C 2-|| In the field range H < H C 2 
at T = we find from Eqs.(l)-(2) 
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Here k„ is the normal state thermal conductivity, 
the effective mean free path l e g = v/(2a(oj = 0)), 
and /it = 4y / 7r(Z e ff/A)(5 2 . Near H c2 (i cx A 2 cx 
H C 2 — H, so that k(H)/k„ m 1 — /i/2 decreases lin- 
early with field. The range over which the linear 
dependence holds, fi -C 1, is smaller the cleaner 
the sample since for ! > ^ we may have /i > 1 
even when S <C 1 . Consequently the linear behav- 
ior down to 0.7H C 2 || indicates a relatively dirty 
sample, in agreement with the estimate of Ref . |3j 
which yielded / ~ 350A, or 1/^q ~ 7. From this 
estimate we obtain T ~ 0.4meV; k xx determined 
with these parameters is shown in Fig.l. 

At low T, for fields below 0.5H C 2, n xx (H) is ob- 
served to decrease with increasing field. || Note 
that the residual linear term in k xx (H = 0,T) is 
only observable if T < <t(cj = 0), which is expo- 
nentially small, in Aq/T, in the Born limit. The 
results of Ref. || therefore suggest that the phase 
shift is close to n/2. But, in the unitarity limit 
cr(w = 0) w V^Aa, and for T ~ 0.4meV at low T 
k xx (H) increases rather than decreases over al- 
most the entire field range. If inelastic scatter- 
ing is strong, T at low T and H may be much 
smaller than its normal state value, in which case 
k xx does decrease but in a narrow field range, see 
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Figure 1. Field dependence of k xx (H,T = 0) for 
T = 0.15A , a = 1.3, and H c2 = 5T. Inset: T = 
0.03Ao, unitarity limit, koo = tt 2 TNqv 2 /6Ao is 
the universal conductivity. Jl[] 



Fig. 1, in disagreement with experiment, suggest- 
ing that the low-field decrease is in the phonon 
contribution to n xx . 
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